Static and Dynamic Properties of 
Inhomogeneous Elastic Media on Disordered Substrate 
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The pinning of an inhomogeneous elastic medium by a disordered substrate is studied analytically 
and numerically. The static and dynamic properties of a D-dimensional system are shown to be 
equivalent to those of the well known problem of a D-dimensional random manifold embedded in 
(D -\- _D)- dimensions. The analogy is found to be very robust, applicable to a wide range of elastic 
media, including those which are amorphous or nearly-periodic, with local or nonlocal elasticity. 
Also demonstrated explicitly is the equivalence between the dynamic depinning transition obtained 
at a constant driving force, and the self-organized, near-critical behavior obtained by a (small) 
constant velocity drive. 



PACS numbers; 46.30.Pa, 81.40.Pq, 64.60.Ht, 
I. INTRODUCTION 

The pinning of elastic continuum in random poten- 
tial has. been a subject of numerous studies in the past 
decadetl. It is related to various phenomena of techno- 
logical importance, while being also of fundamental in- 
terest to the statistical mechanics of disordered systems. 
Considerable efforts have been devoted to understand- 
ing how quenched impurities, influence the transport of 
charge-density wave (CD.W)a"Q, pinning of flux lines in 
type-JL superconductorsB, roughness o£ crystalline sur- 
facesHLl, propagation of invasion frontsQ, etc. Many as- 
pects of these systems can be described either by the 
model of randomly-pinned CDW or by the model of 
randomly-pinned directed manifold ( "random manifold" 
for brevity) . The static properties of the low-temperature 
glass phases of these systems have been studied by a va- 
riety of analytical methods, including renormalizai-iQii.- 
group (RG) analysiaatlpipplica variational methodUO, 
and functional RGEj^tSEJ. Combined with the exact 
ground state structures obtaiped— numerically via effi- 
cient (polynomial) algorithmalaO, these systems are 
perhaps the best characterized glassy system to date. 
Much progress have also been made in understanding 
the nonequilibrium driven dynamics of these systemsEJ. 
In the extreme nonequilibrium limit where thermal fluc- 
tuations can be neglected, it is known that the driving 
force exceeding a critical value is necessary to depin the 
system. A continuous dynamic phase transition occurs 
at the depinning threshold, where the dynamics exhibits 
complex stick-slip motion with "avalanches" of all sizes. 
Such complex dynamics results from the intricate balance 
of elasticity and random pinning forces near the onset of 
motion. They have been chacactecized in great-details 
by a combinjarticm of analyticQO^O, numericE3"Ej, and 
experimenta]E3~E3 methods. 

Another class of pinning phenomena which have at- 
tracted much attention is the tribology of sliding elastic 
bodies interacting via a contact surfacej-This is exempli- 
fied by the Burridge-Knopoff modeEa~E3 describing the 



dynamics along an individual earthquake fault. Other 
examples inphidp boundary layer lubricationES, and stick- 
slip motionc3o of elastic continuum over sticky sub- 
strates. Complex spatio-temporal dynamics are found to 
occur in these sy|Stems also, when they are driven slowly. 
Previous studiesc3 of this class of systems have focused 
on the role of chaos generated by nonlinear dynamics, 
e.g., a velocity- weakening friction. Complexity in these 
systems are not as well understood theoretically. 

Another approach to studying the spatio-temporal 
complexity generated in tribology-like problems is by 
"sandpile"-like autoitiaton models, as pioneered by Bak 
and coUaboratorsO'EJ. The connection between sandpile- 
like models and dynamic critical phenomena has been 
proposed early on by Tang and BakCJ. (In fact, the sand- 
pile modelo itself was motivated from the study of the 
randomly-pinned CDWcZI. Equivalence between many 
aspects of the sandpile model and those of the CDW at 
the depinning threshold has since-ieen verified numer- 
ically by Narayan and Middletonta.) Analogy between 
the geometrical structures of individual earthquake fault 
zones and those ofJiie equilibrium randonp-manifold has 
also been explorecO. In a previous letterEI, we demon- 
strated analytically and numerically the equivalence be- 
tween a certain class of disorder-dominated tribology 
problem and the depinning dynamics of the random mafc 
ifold. The same approach has been used by Fisher et ati 
as a starting point in analyzing earthquake fault dynam- 
ics. Connection between a related interface depinning 
problem and certain earthquake mcjdels have also been 
discussed by Paczuski and BoettcherH In this paper, we 
examine in detail the relation between the tribology-like 
systems and the randomly-pinned CDW/manifolds. Our 
major result, that both the static and dynamic properties 
of a D-dimensional inhomogeneous elastic body embed- 
ded in a D-dimensional random medium are equivalent 
to those of a homogeneous D-dimensional directed man- 
ifold embedded in {D -\- D)- dimensional random medium, 
is found to be very robust, applicable to a wide range 
of elastic media, including those which are amorphous, 
nearly-periodic, with local or nonlocal elasticity. These 



1 



results are relevant to a number of apparently unrelated 
problems, including the enhanced pinning of entangled 
flux Imes, the nonequilibrium freezing of napving vortex 
arraytj, the reptation p£, heteropolymersEj, and align- 
ment of DNA sequenccf£3. We also demonstrate explic- 
itly the connection between the dynamic critical phenom- 
ena obtained from a constant driving force at the dcpin- 
ning threshold, and the self-organized, nearly-critical be- 
havior obtained by a (small) constant velocity drive. The 
latter is an example of "extremal dynamics^ by which 
many self-organized critical phenomena ariscEj. 

This paper is organized as follows. We first introduce 
in Sec. II a discrete D-dimensional model of inhomo- 
geneous elastic system pinned in a disordered medium, 
and derive the appropriate continuum Hamiltonian. The 
statics of the continuum system is examined analytically 
in Sec. III. A. We show that while a perfectly periodic 
system is equivalent to a randomly-pinned CDW system, 
any quenched-in inhomogeneity changes the universality 
class to that of the random manifold. This includes some- 
what surprisingly even nearly-periodic systems with a 
small concentration of quenched-in interstitials and va- 
cancies, or with only quenched-in phonon modes. This 
finding is demonstrated numerically by detailed analysis 
of the model in D = 1. All universal quantities exam- 
ined, including amplitude ratios, are found to be indis- 
tinguishable from those of the one-dimensional random 
manifold, i.e., the randomly pinned directed path in 1 + 1 
dimensions. The driven dynamics of the elastic system 
(with a constant force) is described next in Sec. IV. The 
random manifold analogy is extended to include critical 
depinning dynamics, and demonstrated explicitly for the 
case D — 1. In Sec. V, we study the effect of nonlo- 
cal elasticity mediated by the elastic body in the bulk 
not in contact with the disordered substrate. The dy- 
namics is shown to be analogous to the appropriate ran- 
dom manifold problem with nonlocal elasticity. Finally, 
we compare the dynamics obtained with constant-force 
drive and that with constant-velocity drive. We find the 
two to be nearly equivalent close to the depinning thresh- 
old where the average motion is slow. In order for this 
manuscript to be self-contained, we provide brief reviews 
of the basic properties concerning the statics and dy- 
namics of the CDW/random manifold systems within the 
text. More detailed explanations of the coarse-graining 
procedure and a review of bulk-mediated elasticity are 
provided in Appendix A and B. 

II. THE DISCRETE AND CONTINUUM MODELS 

In this section, we considered a D-dimensional in- 
homogeneous elastic medium, e.g., an amorphous solid, 
which is completely immersed in a Z)-dimensional dis- 
ordered substrate. An example for the D = 2 case is 
a sheet of latex membrane in contact with glassciJ, or 
a randomly-polymerized membrane adsorbed on a sub- 



strate (see Fig. |^). This situation also arises in the pin- 
ning of a rigid vortex array (with quenched-in defects) in 
a thin- film superconductor. Moreover, D = 3 may de- 
scribe the pinning of an entangled vortex line network in 
a bulk superconductor (see below) , and I? = 1 is relevant 
to the jfcptation of a heteropolymer in a disordered gel 
matrixa. 




FIG. 1. A randomly-tethered elastic membrane in contact 
with a disordered substrate. 

We use an irregular array of beads tethered by har- 
monic springs to model the amorphous solid. Let the 
average inter-bead distance be a and the equilibrium po- 
sition of a bead labeled by n be given by R„ G . The 
neighboring beads are connected by harmonic springs 
with appropriately chosen spring lengths (of the order 
a) such that the configuration {Rn} is the unfrustrated 
ground state of the tethered system in the absence of any 
external forces. 

To describe the large scale properties of such an elastic 
medium analytically, it will be useful to adopt a contin- 
uum description. Let the density field of the unperturbed 
system be 

Po(r)-5]<5^(r-R„). (1) 

{"} 

The quenched-in density variation is 5po{r) — Po(r) — p, 
where p = is the average density; it is characterized 
by the correlation function 

Spo{r)Spo{r')=C„{r-r'), (2) 

or the structure factor 5*0 (k) = J d^r Co(r)e''"". Here the 
overbar denotes average over the ensembles of quenched 
bead positions {Rn}. In the case of a single large sys- 
tem, the overbar can be taken as the spatial average over 
smaller subsystems. 




FIG. 2. The equilibnum position (circle) of a bead n is 
given R„; the actual position (black dot) r„ is given by the 
displacement vector u(R„). Note that the set {Rn} is not 
necessarily ordered periodically. 
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Elastic distortion of the discrete system is described 
by a displacement vector u(R„), which denotes the dis- 
placement of the n"^ bead from its equilibrium position 
R„; see Fig. ^. Then the actual position of the bead is 
r„ — R„ + u(R„), and the corresponding density field is 

p(r;u) = ^,5^(r-R„-u(R„)) 

{n} 

V.u)5]5^(r-R„-u(r)), (3) 

{»} 

where we assumed that variations in the displacement is 
small, i.e., V ■ u ^ 1. Comparing Eq. (||) with Eq. (|^), 
we see that 



p(r, u) = (1 - V • u) • po(r - u) 

« p(l V • u) + Spoir - u) 



(4) 



to leading order in (V ■ u). In the continuum description, 
details of the elastic medium are contained completely in 
the term 5po{r), through the correlation function Co(r). 

In the absence of any external forces, the energy of the 
tethered system is invariant upon a constant displace- 
ment u — > u+const. For a statistically isotropic medium. 



the energy of small e 
by the classical form 



stic distortion is then simply given 



Ho 



V-u) 



C66 



(V X U)2} 



(5) 



where Cn and Cgg are the bulk and shear elastic modulill. 

Consider now the situation where the elastic medium 
is in contact with a disordered substrate, modeled by a 
Gaussian random potential V(r) with zero mean and a 
variance 



[y(r)l/(r')] = Ay^^(r-r'). 



(6) 



Here, [...] denotes average over the ensemble of sub- 
strates. The interaction of the medium with the sub- 
strate is described by a pinning energy. 



H 



pm 



where 



F(r„) = y"d^rp(r,u)F(r) 
y"d^r{ -p(V-u)t/(r) + W^(r,u(r))}, (7) 



E 

{"} 



W{r,u)=6po{r-u)V{r) 



(8) 



and the term pV{r) is neglected as it produces merely an 
overall energy shift. Note that W{r, u) depends explicitly 
on u where as other terms depend only on Vu. Thus, 
only W breaks the translational symmetry u u+const, 
and is responsible for providing the pinning phenomenon. 
The interacting system, characterized by the effective 
Hamiltonian 

n^Ho + J d^r{-p(V-u)l^(r)+I^(r,u)}, (9) 
will be analyzed in detail in the next sectionS. 



III. THERMODYNAMIC PROPERTIES 
A. Theoretical Considerations 

1. The Periodic Medium 

If the tethered system does not contain any quenched- 
in defects, then the intrinsic density variation is periodic, 
i.e.. 



Spo{x) p ^ cos (Ki • x) , 



(10) 



where K^'s are the reciprocal lattice vectors. Eq. (||) in 

this case reads 

HcDW = ^0 + y rf'^rj - p (V • u)F(r) 

+ ^py(r)cos[K,.(r-u(r)]}. (11) 



The model (^T|) is the ZJ-component generalization of 
the randonilv-pinned charge- density wave (CDW) in D- 
dimensionsE^til. The equilibriurri|-jp£nperties of this class 
of systems have been well studiedBliy'tHl: It is known that 
the disorder is irrelevant in Z? — 1, where the system 
behaves like a Gaussian chain. In D >^2_the system 
( |rT| ) is glassy at any finite temperaturesEdO. The glass 
phase is described by two critical exponents: the ther- 
mal exponent 9 characterizing typical fluctuations of the 
free energy landscape AJ^ ~ for elastic distortion 
over the length scales L ^ a, and the "wandering" ex- 
ponent C characterizing fluctuations of the displacement 
fleld, Au ~ L^, in the low free energy state(s). The 
results 6 = D ~ 2 and C, = O(log) is believed to be ex- 
act for 2 < D < 4. Right |ijT_|Q-.= 2, the situation is 
somewhat more complicatedcl'EjS. The disorder is ir- 



relevant for T > Ta 



value of the 
scopic modell 
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The precise 



I temperature T„ depends on the micro- 



Below T, 



9' 



dimensional system 
is in a "marginal" glass phaseHEH'Ell characterized by 
logarithmic anomalies in AT and Am. 



2. Strongly-disordered Medium 

If the tethered system contains a finite concentration 
of quenched-in dislocations and/or disclinations, then the 
density variation 5po is no longer described by (|l^), and 
one must characterize 5po statistically through the corre- 
lation function Co(r) or the structure factor 5*0 (k). Con- 
sider the simpler case where 5*0 (k) is liquid-like, contain- 
ing no Bragg peaks at finite k. This may describe, for 
example, a completely disordered film such as a rubber 
sheet, or a randomly-polymerized liquid membrane. The 
model defined by Eqs. (^ and now contains two kinds 
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of (mutually uncorrelated) disorders in the effective ran- 
dom potential W{r,u.), making it difficult to solve sys- 
tematically. In fact, straightforward application of the 
replica trick to the model (^ immediately leads to difh- 
culties as one has to integrate out the disorders twice. 

However, if we simply treat W{r, u) as an effective dis- 
order and examine its moments using (^) and (|^) , we find 
it has zero mean and a variance given by 



revealed by considering higher moments of W, for exam- 
ple, W^2(r, u) = W^{r, u) - [W^]. We find 



[W(r, u)T4^(r', u')] = AvS"{r - r') Co(u - u') 

« A(5^(r-r')(5^(u-u'), (12) 

with A = Ay 5*0 (k 0) for short-range correlated func- 
tion Co(r). 

Note that the form of the correlator (|2|) is the same 
as that of an uncorrelated random potential in the 
space (r, u) G 5R^'*"'°. It is then tempting to interpret 
the system (^ as a ZJ-component, D-dimensional "di- 
rected manifold" u(r) embedded in an effective {D + D)- 
dimensional random potential M^(r, u). The latter is 
an example of the so-called "random manifold" prob- 
lem which is encountered in a j uid j e . ya riety of context 
involving quenched randomnessB'PrPirT. Properties of 
the random manifold (RM) have been characterized in 
detail and are briefly summarized here: It is known that 
an d-component RM in £)-dimension is asymptotically 
described by a glass phase at any finite temperatures if 



2D>d{2- D) 



(13) 



Outside this regime, the glass phase is still obtained-if the 
temperature is below a certain critical temperaturetZI T > 
0. Like the randomly pinned CDW, the behavior of the 
RM in the glass phase can be described by the thermal 
exponent 9 and wandering exponent ^. There is an exact 
exponent identity 9 — D — 2 + 2( linking the two expo- 
nents for all d andpD. The exponents arc known exactly 
for the special casal3 D — I and d — 1, with ( = 2/3 and 
9 = 1/3. There are also strong bounds on the exponents, 
with (4-D)/4 < C < max{(4-i:i)/(4-f d),(2-L>)/2}. 
[Note that C ^ as _D — > 4^, reflecting the fact that 
D = 4 is the upper critical dimension of the problem.] 
Numerically, the exponents foij-izjarjous d and D have been 
determined to good accuracylljilZl. The results are ap- 
proximately summarized by the expression 



c 



2(4 -D) 
8 + d 



(14) 



which was motivated, by a functional renormalization- 
group considerationE3. 

What does the random manifold problem have to do 
with the problem at hand ? Even though the sec- 
ond moment of W{r, u) is of the same form as that 
oi a {D + _D)-dimensional random potential, W itself, 
being the product of two D-dimensional random func- 
tions, certainly cannot be truly be a short-range cor- 
related {D + £))-dimensional random potential. There 
must be long-range correlations, the forms of which are 



[1^2 (r, u)I¥2(0, 0)] = A6'{r)S^{u) + 2(5^ (r) + 2(5^ (r - u), 

(15) 

indicating correlations along the "directions" of constant 
r and constant r — u. Such correlations are of course not 
surprising given the form of W^(r, u) in Eq. (^). Thus, 
a more accurate model of the effective random potential 
should include a superposition of short-range correlated 
and long-range correlated random potentials, e.g.. 



W{r, u) = Wo(r, u) + VKa(r) + Wb{r- u) 



(16) 



with Wo{r,u) being truly a short-ranged {D + D)- 
dimensional potential described by the correlator 
(|T^). The random potentials Wa and Wb generated 
from the higher moments of W have the correlator 
[V7,(x)W^,(x')]~<5„-^(x-x'). 

Clearly, Wa (r) has no effect on the behavior of the sys- 
tem as it does not involve u. Since |u| <C |r| {C, < 1), the 
u-dependence in Wb(j — u) is perturbatively irrelevant, 
and we conclude that Wb does not play a significant role 
either in the limit of weak disorders. Another way to gain 
some intuition of the correlated potential is through an 
example in _D = 1. In Fig. p|, w e illustrate the form of the 
effective random potential (Eq) in Z) = 1, where the "ran- 
dom manifold" is usually referred to as a "directed path" 
(directed horizontally in the direction r of Fig. §), and 
the correlated potentials are like the "columnar defects" 
encountered ipr-|t|lie pinning of flux lines in high- Tc su- 
perconductorg^JlId. These columnar defects are known 
to be strongly relevant if oriented in the direction of the 
directed path (the r-direction in Fig. ||). However, they 
are irrelevant if tilted away by a slope exceeding a finite 
threshold given by the strength of the disorder. The cor- 
related potentials Wa and Wb are just the D-dimensional 
generalization of the "tilted" columnar defects. Their ir- 
relevance can be established more rigorously by extend- 
ing, for example, the analysis of Ref.Ej to D-dimcnsions 
and will not be pursued here. 




FIG. 3. Schematic illustration of the effective random 
potential (^^ in D = 1, for a path directed in the horizontal 
r-direction: The correlated component is shown as two sets of 
"columnar defects", pointing along the directions r = const 
and r — u — const; the uncorrelated component is shown as 
"point defects". 
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For strong disorders, a large distortion with u r is 
admissable by our model defined by Eqs. (^) and (||): 
Such a distortion would be favorable if the elastic en- 
ergy cost (per volume) is more than compensated by the 
disorder energy gained. The latter is given by the or- 
der of variations in Wb and is finite. Thus, our model 
can in principle display a phase transition to a "local- 
ized phase" where u ^ r, similar to what was found 
by Hatano and Nelson in the context of non-Hermitian 
quantum mechanicalj. However, the formation of this 
localized phase would require different "beads" of the 
manifold all to lie with a finite volume of the substrate. 
This is clearly unphysical and is prevented in practice by 
any excluded- volume interaction between the beads. 

Based on the above analysis, we conjecture that the 
long-ranged correlations in W(r, u) are irrelevant for ar- 
bitrary disorder strengths, and the problem of a random 
elastic medium on a disordered substrate belongs to the 
same universality class as that of the random manifold 
with d = ZXjO. Since the condition ( p^ ) is always sat- 
isfied for d = D, we expect our system to be glassy at all 
finite temperatures, with the exponents 



c 



2{A-D) 
8 + D 



and 



0(2 + D) 
8 + D 



(17) 



upon adopting the approximate formula (p^. 

Our conjecture was presenteiior the case D = 1 previ- 
ously in a short communicationtHl. In Sec. III.B below, we 
present results of an extensive numerical study. We find 
that all universal aspects of the D = 1 problem mea- 
sured, including amplitude ratios, agree quantitatively 
with those of the 1-1-1 dimensional directed path in ran- 
dom media, thus vetLfying our conjecture for D = 1. Very 
recentivpZeng et aO applied the min-cut-max-flow algo- 
rithmllJ'Ea to investigate the D = 2 version of the model 
defined by (||) and (||), but with only one component of 
the displacement field u. They found an exponent value 
C, ~ 0.42, which is consistent with that of the random 
manifold with d = 1 and £> = 2 (C = 0.41 ± 0.01), as ex- 
pected according to our conjecture. For a 2-dimensional 
random elastic medium on a 2-dimensional substrate, we 
predict that C « 0.4 and 6 w 0.8 as for a 2-component, 2- 
dimensional directed manifold in 4-dimensional random 
medium. An entangled flux Jipe array, e.g., a "poly- 
mer glass", may be modeledLJ by a 2-component, 3- 
dimcnsional random elastic system for time scales up to 
the distanglement time r ~ e^^/"^, where Ux is the flux 
cutting energy estimated to be of the order of 10 times the 
melting temperaturdl3. The corresponding exponents are 
thus expected to be C « 0.2 and 9 « 1.4 in the elastic 
regime. 



This could be the case of a rigid vortex array on a thin 
film superconductor ,|-t|he defects being frozen- in vacan- 
cies and interstitialsE^I. Since a low density of vacan- 
cies/interstitials does not destroy the crystallinity of the 
elastic system, Bragg peaks in 5*0 (k) still leads to a CDW- 
type interaction ( |l l| ) when the elastic medium is placed 
in contact with the disordered substrate. However, the 
presence of quenched-in vacancies and interstitials also 
gives rise to large-scale density variations, manifested by 
a nonzero component of 5*0 (k — > 0), which leads to the 
effective [D + L')-dimensional random potential as de- 
scribed by the correlator ([l^). Thus, the nearly-periodic 
elastic system is subject to both the CDW-type and the 
RM-type disorders. What is the outcome of competition 
between these two kinds of interactions ? A simple power 
counting reinforces the intuitively obvious result that the 
RM interaction is relevant in the CDW phase, while the 
CDW interaction is not relevant in the RM phase. We 
thereby conclude that even the nearly-periodic elastic 
system belongs to the RM universality class. 

Our analysis indicates that any finite concentration 
of quenched-in defects is sufficient to change the CDW- 
type pinning of a perfectly periodic system to the much 
stronger pinning of the random-manifold universality 
class. This includes a low concentration of quenched-in 
interstitials and vacancies, which by themselves do not 
destroy the crystallinity of the elastic medium. Such a 
result may be rather surprising at a first glance, since 
the existence of periodicity of an elastic medium is usu- 
ally associated with the CDW universality class. But 
this notion is incorrect. What is responsible for the 
CDW universality class in the pinning of an usual pe- 
riodic solid is a relabeling symmetry of the underlying 
discrete system, e.g., the energy of a configuration of 
beads described by {Rn} is the same as that described 
by {Rn+const} (up ta-boundary effects). This relabeling 
symmetry is brokenEII given a finite concentration of in- 
terstitials/ vacancies, since different beads are no longer 
equivalent. Thus, asymptotically, the system is con- 
trolled by the RM fixed point. 

Lastly, we mention that even if the medium contains no 
topological defects at all (not even interstitials and vacan- 
cies), some quenched-in "phonon modes" may already be 
sufficient to induce the RM behavior: Let the equilibrium 
positions of the beads labeled by (ni,n2) be Rni.na = 



+ w„ 



where TZn 



denotes the points of a 



periodic lattice, and w is the quenched-in phonon modes 
characterized by [(w(x) — w(x'))^] = A(x — x').|-Then 
the correlator of the effective pinning energy W iala 



[M^(r,u)W^(0,0)] = Ay 5' 



'(r)^cos(K, 



u e 



2K^ 



3. Nearly-periodic Medium 

We next turn to the case of a nearly-periodic, tethered 
system with a low concentration of quenched-in defects. 



where 's are again the reciprocal lattice vectors. Since 
\W{r, \i)W{Q, 0)] is short-ranged correlated in u as long 
as A(u) diverges for large |u| (including logarithmic di- 
vergence) , we see that the RM universality class is recov- 
ered even for quenched-in phonon fluctuations in Z? < 2. 
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B. Transfer Matrix Studies 

We test our predictions by performing numerical stud- 
ies of a one-dimensional bead-spring system correspond- 
ing to the D = 1 case of the randomly-tethered elastic 
medium considered above. The one dimensional system 
is chosen since its thermodynamic properties can be ob- 
tained in polynomial (iV^) time using the transfer ma- 
trix methodlij, and also the thermodynamics of the cor- 
responding (l-l-l)-dimensional problem of a diji^cted path 
in random media (DPRM) is known exactlyll3. Thus a 
quantitative comparison can be made. 




FIG. 4. Longitudinal displacement variation of the beads 
of the random chain in (a) is represented as a directed path on 
a (1 -I- l)-dimensional lattice in (b). The equilibrium spring 
lengths are ai — 4, a2 = 5, and aa — 3. The equilibrium 
positions R„ — are indicated in (a) as solid lines. 

The open circles indicate the allowed bead positions, given 
the SOS-like restriction on u„. 



We consider the following discretized one-dimensional 
problem: A chain of + 1 beads (labeled sequentially 
by n G [0, iV]) is placed on a one-dimensional lattice of 
unit lattice spacing. Each bead n (except for n = 0) is 
connected to its nearest neighbor n — 1 by a harmonic 
spring. All springs have the same spring constant 7, but 
the equilibrium length a„ is an integer drawn randomly 
from the interval [5, 15] (in units of the lattice spacing), 
such that the mean spring length is a = 10. The equilib- 
rium positions of the n*^ bead is i?„ = J2m=i '^"i if we 
fix the n = end of the chain at the origin. To speed 
up the numerics, we apply an SOS-like restriction and 
allow the springs to be compressed or stretched by at 
most one lattice unit. (We have verified that allowing for 
large excitations does not affect scaling properties of a 
long chain.) Then the chain configuration is given by the 
position of each bead n 



(18) 



where m„ e {0, ±1,±2,... ,±n} is the "displacement 
field" , with the restriction (u„+i — m„) e {0, ±1}, and the 
energy of each spring is ^{un+i — — {O, }■ Since 
each chain configuration is uniquely specified by the set 
of numbers {un\ = (mi,U2, . . . ,un), we may represent 
the chain as a directed path in 1 + 1-dimcnsions, with 
"transverse" coordinate {u„}. The mapping is illustrated 



in Fig. ^. The allowed bead positions are marked by 
open circles in Fig. ^(a) and shifted upwards for clarity. 
The chain with -t- 1 beads has a total of 3^ different 
configurations. One of these configuration with {it„} = 
(—1,0,1) (as shown in Fig. ^a)) is represented by the 
full line in the directed path representation of Fig. 0(b). 



(b) uniform 
springs 





FIG. 5. Gray-scale plots of potential Wn(u) — V{R„ +u), 
Rn = im, for a chain with A'^ = 100 elements, (a) Ran- 

dom springs with a„ uniformly distributed in interval [5, 15], 
and (b) uniform springs with a„ — 10. Lighter shades corre- 
spond to larger values of W. 



(a) random 
springs 




FIG. 6. Collection of the locally optimal paths for disorder 
potential W shown in Figs, ^(a) and (b). The bold lines 
represent the global optimal paths. 

Next, we specify the energetics of the system. Dis- 
order of the substrate with which the chain is in con- 
tact is given by the function V{r), drawn independently 
from a Gaussian distribution of zero mean and unit vari- 
ance for each lattice point r. The spring constant used is 
7 = 0.2 such that the typical "spring energy" is 0.1, much 
smaller than the disordered "potential energy" ~ 0(1). 
This choice of parameters facilitates a quick approach to 
the asymptotic (glassy) regime. In term of the displace- 
ment variable u„, the random potential can be written 
as Wn{un) = V{Rn + u„), which depends explicitly on 
the two sources of randomness, V{r) and It is illus- 
trative to plot Wn{u) in the two dimensional space (n, u) 
(Fig. |^(a)). For comparison, the case of uniform springs 
(with an = 10 for all n's) is shown in Fig. H(b). We see 
that the periodic feature of W shown in Fig. @(b) is ran- 
domized by the random springs. However, correlations 
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in a slanted direction is still clearly visible in Fig. ||(a) 
(cf. Fig. I). 

Effect of the different potential W on the directed path 
can be illustrated by differences in the morphology of the 
"local optimal .paths" , i.e., the collection of the lowest en- 
ergy pathst^ilj of length N, connecting the starting point 
at the origin and all possible ending points at {N, u) . In 
Figs, ^(a) and (b), we plot the optimal paths for the 
realizations of W shown in Figs. ||(a) and (b) respec- 
tively. The bold line corresponds to the global optimal 
path, the lowest energy path among all the local optimal 
paths. Note that the local optimal paths of the system 
with uniform springs (Fig. ^(b)) are very regularly ar- 
ranged, with the distance between neighboring branches 
almost constant and approximately equal to the equi- 
librium spring length a. On the other hand, the local 
optimal paths for the random chain are arranged much 
more irreg ularly (Fig. |(a)). There are for example large 
islands of high energy regions which the paths avoid, sim- 
ilar to what was found previously foTj^he directed path 
in 2-dimensional random potentiaO'tj. Thus, we can 
view these paths configurations as an indication of pos- 
sible equivalence between the statistical properties of a 
random chain on one-dimensional random substrate and 
a directed path in 2-dimensional random medium. 

The statistics of the optimal paths have been investi- 
gated in Ref. |5^, where we presented numerical results 
obtained from the zero-temperature transfer matrix so- 
lution of systems with N — 4096, averaged over 2000 
independent realizations of V{r) and {a„} pairs. Fluc- 
tuations in the end-to-end displacement it^ and the to- 
tal energy E*{N) of the global optimal path {u*} indi- 

cate scaling behavior with [(u^)^] ^ [u^] ~ N'^^^ and 

[{Elf)'^] - [E"^] 7V^/^, both of which are charactecisr 
tic of the 1-1-1 dimensional DPRM universality classll3. 
Here, we describe extension of the zero-temperature cal- 
culations to finite temperatures. 
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FIG. 7. Log-log plot of the sample-to-sample free en- 
ergy fluctuation (AJ^)'^ (N) for various chain lengths A''. Thin 
straight line shows the predicted asymptotic scaling. 




FIG. 8. Log-log plot of the end-to-end displacement fluc- 
tuations for different chain length A^, at the temperature 
T = 1. Also shown is the connected thermal fluctuation Su%. 



The Boltzmann weight WnIu) for paths connecting 
their origin at (0, 0) to their end |P'f|^ts at {N, u) can be 
obtained recursively according tolj'Ea: 

Wn+i(u) = co{n,u)Wn{u) 

+ci{n,u){Wn{u-l) + Wn{u+l)Y (19) 

where co{n,u) — e^"^"^/^, ci{n,u) = co{n,u) e^''^^^'^\ 
with the "initial condition" Wo{u) ~ <^n.o- The partition 
function Zn is obtained as Zn — W„(u), and the free 
energy is J-n = —TlnZn- The thermal averaging of an 
observable On {u) is given by 



(20) 



The sample-to-sample free energy fluctuation (AjT)^ 

-2 



[J-jf] — [J'n] is shown in Fig. |^. The results are obtained 
at temperature T = 1.0, from systems with A^ = 1024 
averaged over 4000 different realizations of random V{r) 
and {an}- It is seen that (AjF)^(A^) approaches the 
asymptotic scaling form N'^/'^. We also computed the 
finite-temperature end-to-end displacement fluctuations 
[(u^)] and [{un)'^]- As shown in Fig. |[ the approach 
to the expected N^/^ behavior is clear. More inter- 
estingly, note that the disorder average of the ther- 
mal fluctuation itself, i.e., the connected second moment 
Suj^ = [{ulf) — {un)"^] in fact scales as A^^ as if random- 
ness is not present. This behavior is esnectad of the 
DPRM due to a statistical tilt symmetryEJ'oEl. Such a 
symmetry is obviously not present in the bare potential 
Wn{u) shown in Fig. ||(a). The scaling law on found 
is therefore another evidence indicating the irrelevancy 
of the slanted correlation in Wn{u) and the dominance 
of the I+l dimensional DPRM behavior. 
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FIG. 9. Log-log plot of gN{T) as a function of temper- 
ature, averaged over 10'' different realizations of disorder. 
The high temperature behavior approaches which is 

expected of the 1-1-1 dimensional DPRM. 
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FIG. 10. Probability distribution for the free energy, plot- 
ted for the rescaled variable, cJjT = (^J^ — [T]^ /AT{N), at 
T = 1.0 for system sizes ranging from TV = 32 to = 1024. 
Data were collected from 4096 samples. 



similar to that of the (1-1- l)-dimensional DPRM which 
has no finite temperature phase transition^. Wc there- 
fore conclude that in the thermodynamic limit N ^ (x, 
there is no finite temperature phase transition, and the 
large scale behavior of the random chain is always de- 
scribed by the zero-temperature (glassy) system. 

We have also calculated the full free energy probability 
distribution P for different lengths N. Figure |l^ shows 
data collapse of Pn{5^) for = 32 to = 1024, col- 
lected from 4096 samples, with 5T = {F - /AT{N) 
being the dimensionless measure of free energy variation. 
The distribution is asymmetric. The skewness 73 and 
kurtosis 74 of this distribution are plotted in Fig. [Tl](a) 
and Fig. |l^(b) respectively for different TV's. Both 73 
and 74 are found to approach the corresponding universal 
numbeips,(dashed lines) known for the (H-l)-dimensional 
DPRME3: 73 = -0.296 ± 0.028, 74 = 3.16. 

Putting together all the numerical results presented 
in this section, we see strong evidence supporting our 
expectation that the thermodynamics of a random chain 
on the disordered substrate and the (l-l-l)-dimensional 
DPRM indeed belong to the same universality class. 
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Our results suggest that the finite temperature behav- 
ior of the system probed is dominated by the T = 
fixed point. This result is not immediately generalizable 
to all T. Naively, one might think that at sufficiently 
high T such that the thermal fluctuation Sun exceeds 
the quenched variation in _Rjv, then the effect of random 
spring length may be washed out. To investigate the pos- 
sible existence of a finite temperature jphase transition, 
we compute the dimensionless quantityCj 



3.50 



9n{T) = [{unY^ 



IK)] 



(21) 



which vanishes for T 00 while g — » 0(1) for T ^ 0. 
If there is a phase transition at some finite temperature 
Tc, the expected finite-size scaling in the critical region 
would be gAr(T) = g{N^I''{T - Tc)), where z/ is the cor- 
relation length exponent. Therefore, the curves of g{T) 
for different A^'s should all intersect each other at Tc if a 
finite temperature transition exists. The numerical data 
for (/Ar(T) in the temperature range T — 0.1 to 3.0, for 
system sizes from = 4 to iV = 128 are shown in Fig. ||. 
No indication of curve crossing at T > is found. In fact, 
the high temperature behavior of gN{T) approaches l/T^ 
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FIG. 11. Skewness 73 and kurtosis 74 of the free energy 
distribution for different chain lengths A*'. Dashed lines are 
the corresponding values for the (1 -|- l)-dimensional DPRM: 
73 « -0.296 and 74 ~ 3.16. 



IV. DRIVEN DYNAMICS 



A. Theoretical Considerations 
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1. Driven CDW and Random Manifolds 

The zero-temperature driven dynamics of the tethered 
system is of interest to the study of tribology and to 
understanding the nonequihbrium dynamics of vortices. 
From the correspondence between the thermodynamic 
properties of the randomly-tethered elastic system and 
the CDW/RM systems, it is tempting to conjecture that 
the correspondence persists also for the dynamical prop- 
erties. Before we provide evidences in support of this 
generalization, let us first review the known dynamics of 
the driven CDW/RM systems. 

Consider first the simplest zero-temperature driven dy- 
namics of the irajBdomly-pinncd (one-component) CDW 
in Z?-dimensionl3'u, given by the equation of motion 

fio^dtu{r, t) = -—ncDwM + F (22) 

where /zq is a bare frictional coefficient, and F is the 
driving force. For F below some threshold force Fc, the 
average velocity v = {dtu) is zero. [Here, (...) denotes 
temporal and spatial average.] Upon approaching the 
threshold from below, the dynamics (e.g., response to 
perturbation) becomes very "jerky". It is consisted of 
a series of "avalanchesll, whose (linear) size i obeys a 
power-law distributionc3, 

Pr{i> s) ^ s-^pis/O- (23) 

In Eq. (|2^), ^ is the correlation length of the system, 
p{x) is a scaling function which is constant for x 1 
and drops off sharply for a; 3> 1. The correlation length 
diverges as ^ '--^ {Fc — F)~'^ as F F^ . The mo- 
tion becomes continuous for F > Fc due to overlapping 
avalanches. There the interface acquires a finite veloc- 
ity with V ^ {F — FcY^ , similar to the emergence of the 
order parameter in a critical phenomenon. These expo- 
nents have been computed|-|by a functional rcnormaliza- 
tion group (FRG) analysisQ to first order in e = 4 — D, 
with = 1/2 and /3 = 1 - e/6. The one-loop FRG 
results are found to be consistent with extensive nu- 
merical simulations, of the driven CDW in various spa^. 
tial dimensionsES^O. For example, Myers and Sethnacj 
find the one-dimensional driven CDW (e ~ 3) to have 
1/ « 0.4 ± 0.1 and [3 « 0.45 ± 0.05.^, 

Similar depinning phenomenon^ is obtained for the 
zero-temperature driven dynamics of the d-component 
random manifold u(r), whose equation of motion is 

Aio'atu(r,i) = -— 7^rm{u} + F, (24) 

where Hrm refers the random manifold Hamiltonian (|^) 
with a truly random {d -\- iD)-dimensional potential W , 
and F is again the driving force. A continuous depinning 
transition occurs at a critical force |F| = F^ where the 
system exhibits avalanches with power-law distribution 
as in (|2^). In the vicinity of the critical point, the ex- 
ponents V and (3 can be defined as before. [Here, (3 is 



the exponent describing the onset of the parallel com- 
ponent of the velocity, v = (9tu) • F.] As the driving 
forcCpi' breaks the isotropy of the system, it is conve- 
nientcil to divide the displacement u into components 
parallel and perpendicular to F, with U[| = u • F and 
= u — U|| F. The depinning transition can then be 
understood in terms of the critical fluctuations in M|| and 
, given by the correlation functions 

((«ll(r,t)-«ll(0,G))') =|r|2x,,g||(|t|/|r|^«) (25) 
((u^(r,t)-u^(0,G))') ^\v\^^-g^{\t\/\vY-) (26) 

at |F| = Fc- They are characterized by their respec- 
tive roughness exponent X||.± and dynamic exponent 
^11, ± I from which all other exponents can be obtainedEJ. 
For example, v= 1/(2 - x\\), P = v {z\\ - x\\) ami 
K = Z) — 2 + X||"- It was shown by Ertas and Kardar£il 
that the scaling properties of M|| arc the same as those of 
the one-component system (RM—ssath d = \) which has 
been solved by the FRG methocHEl, with X|| = e/3 and 
Z|| = 2 — 2e/9 to_Sjst order in e = 4 — Z?. Ertas and 
Kardar also foundEil X-L — X\\~ % ^i^d zj_ — zy + v^^. 

Numerically, one finds-for the one-dimensional inter- 
face in (l-|-l)-dimensiong23 (i.e, d = iD = lore = 3) that 
X « 0.97 ± 0.05, V « 1.05 ± 0.1, and (3 « 0.24 i,0.1. For 
the two-dimensional interface in 3-dimensiona2il (d = 1 
and e = 2), the exponents are x ~ 0.67 ± 0.03 and 
V w 0.75 ± 0.05. These results are all consistent with 
the one-loop FRG predictionfOnJ. In addition, a differ- 
ent exponeiit is found for the finite size scaling of the 
roughncsscZl, with 

{{u{L,t)-u{Q,t)Y)^L^^--, 

with xfs ~ 1-25 in D = 1. 

2. Driven Tethered Network 

We now return to the randomly tethered elastic net- 
work defined in Sec. II, and study its motion in the pres- 
ence of a uniform driving force F parallel to the substrate. 
We start with the deterministic and purely dissipative dy- 
namics of the discrete bead-spring system. The equation 
of motion, in terms of the displacement vector u(R„) for 
the bead n, has the form 

^lo^^tVL{Iin,t) = - ^ ^"'(R„ - R„') • U(R„',0 

{"'} 

-— l/[R„ + u(R„,t)]+F, (27) 
ou 

where the kernel describes the spring forces exerted 
by all the beads {n'} connected with the bead n. 

As we will show in Sec. IV. B, the behavior of this teth- 
ered network is very similar to that of the driven CDW 
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and RM just described: At small driving forces the sys- 
tem is completely pinned by the random potential. As 
the force increases above some threshold value Fc, the 
system starts to move with a nonzero average velocity, 
V ^ {F — Fc)^ ■ The behavior near the depinning tran- 
sition is characterized by a diverging correlation length 
^ ^ a as in usual critical phenomena. 

To study the depinning phenomenon, it is useful to 
coarseorain the equation of motion for the discrete sys- 
tem (^7|). The procedure is described in Appendix A for 
a 2D system. In term of the coarse-grained displacement 
field u, the equation of motion becomes 

fio'dMv, i) = (1 - V • u)f (r, {u}) + fpin(r, u) (28) 

where 

f(r, {u}) = g-\v - r') • u(r',t) - VV^<(r) + F 

(29) 

results from the straightforward coarse graining of the 
right-hand side of Eq. (p7|), with V<(r) being the slowly 
varying part of the substrate potential V{t). Since the 
beads in the tethered systems are connected only to other 
beads in their vicinities, the kernel is local. Its 
Fourier transform reads in component form 

Gy-^ (k) = CQQSijk^ + (cii - cee)hkj, (30) 

where we have again neglected the spatial variations in 
the elastic moduli. Finally, the u-dependent pinning 
force in (|8|) is 

fpin(r,u) = -V^S(r).^po(r-u)/p, (31) 

where Vy{r) are the Fourier modes of ^(r) close to the 
inverse bead spacing. 

For the system with uniform springs, density variation 
is given by (p^). The equations of motion ( psf ) - ( |3l| ) 
are then similar to the ones describing the randomly- 
pinned (Z?-component) CDW. This is expected given the 
thermodynamic properties described in Sec. III. Com- 
pared to the "usual" equation of motion for the driven 
CDW, fi^^dtu = -^SHcDw/Su + F, Eqs. -j|l|) con- 
tain addition terms such as 'VV and (V ■ u) f . These 
terms have .repently been introduced on phenomenologi- 
cal groundsOlp, Here we find that they can be obtained 
systematicall5itZI from a coarse-graining procedure (Ap- 
pendix A). What effects do these additional terms have 
in the vicinity of the depinning transition ? The term 
Vy4s u-independent and hence does not provide pin- 
ningEa. The term proportional to (V • u) is dynamic 
in origin. Kinetically-generated terms of this kind, in- 
cluding other terms with higher powers in (V • u), can 
drastically affect the behavior of the system in the limit 
of strong drive where f is large. However, they are ir- 
relevant iiL±he vicinity of the depinning transition where 

dtu,t^(m. 



The introduction of random springs destroys period- 
icity in Spo and one must describe the random pinning 
force fpin(i", u) statistically through the correlation func- 
tion Co as done in the equilibrium case. We find 

[/pi„..(r,u)/pin,,(0,0)] = ^ (t) '5«^'5W^o(")' 

« ^6{r)d^djCo{u) (32) 

which is short-range correlated in both r and u for short- 
ranged correlated 6pQ. As in the static case analyzed 
in Sec. Ill, higher moments of fpin are long-range corre- 
lated but irrelevant. The critical dynamics with short- 
ranged pinning foijces is then in the universality class 
of the driven RME3. In fact, the pinning force ap- 
pears as if generated directly from the effective poten- 
tial M^(r, u), i.e., fpm « —SW{r,u)/6u. We thus con- 
jecture that the critical depinning dynamics of the D- 
dimensional randomly-tethered elastic network on disor- 
dered substrate is in the same universality class as the D- 
component, £)-dimcnsional ditected manifold in (D + D)- 
dimensional random mediumEd. 



B. Numerical Simulations 

In this section, we report detailed numerical simulation 
of the critical dynamics of the driven one-dimensional 
bead-spring system whose equilibrium properties were 
described in Sec. III.B. The zero-temperature response 
of the random chain to an external driving force is ob- 
tained by a direct numerical integration of the over- 
damped equation of motion, 

dtrn{t) = 7 [rn+i{t) - 2r„(t) + r„_i(i) 

~ian+i-an)]-V'ir„) + F (33) 

where r„ is the position of the n*"^ bead, a„ is the equi- 
librium length of the n*^ spring and V — —dV/dr. 
Expressed in term of the displacement field Un(t) — 
r„(<) - i?„, where i?„ = Sm^iam, Eq. (|3|) is just 
the one-dimensional version of Eq. (g7|). (The bare fric- 
tional coefhcient fig is set to unity here.) We integrate 
Eq. ( ^ ) in the simple Eulerian manner by discretizing 
time. The positions r„'s are kept as continuous variables. 
The random potential V{r) is constructed by a series of 
(quadratic) bumps and valleys, centered on a lattice with 
unit spacing, i.e., at j = {0, ±1, ±2, . . .}. The range of 
each bump/ valley is Rq < 1/2, and the amplitude Vj of 
the bump/ valley centered at site j is drawnly randomly 
from the interval [—1, 1]. Thus, 

^(0 = E y [('^ - ^■)' - ^o] 0(^0 - |r - J I) . (34) 
j 

Since in Eq. ( p3[ ) there is no term which forbids local 
crossings among the chain elements, in the evolution of 



10 



r„(i) we explicitly put additional restriction eliminating 
such moves, so that r„+i > r„ always holds. This speeds 
up the dynamics and does not affect the asymptotic scal- 
ing behavior as we verified. The simulations were run 
on different system sizes with 7 = 0.1, Rq — 1/2, and 
spring lengths {a„} chosen randomly from the interval 
[5,15]. Various time step sizes were used, ranging from 
At = 0.5 to At = 0.01. For all the results reported, we 
always checked that twice smaller At does not lead to 
significant differences. Open chain boundary conditions 
were imposed by introducing two fictitious beads, with 
r-i{t) = ro{t) - flo, and rpf+iit) = rAr(t) + un+i sup- 
plementing Eq. (33). As initial conditions, we take each 



spring to be either compressed or stretched, within 50% 
of its equilibrium length. The chain is then released in 
the random environment described by Eqs. ( ^3| ) and (34). 
Each bead is pulled by a constant force F which is the 
only parameter we vary. 

By applying a strong enough force, the system starts to 
move with a velocity which after some initial fluctuation, 
reaches its stationary value 



N 



(35) 



Here (• ■ •) represents averaging over long time, which is 
very much needed in the vicinity of the depinning tran- 
sition where the chain motion becomes very jerky as we 
illustrate below: In Fig. |l2|, we show the bead trajectories 
rn{t) for a chain with = 16 beads. For clarity, only 
the trajectories of beads with indices n = 0,4,8,12,16 
are shown. The driving force is F = 0.2 in Fig. p^(a). 
After some initial movements, all the trajectories become 
independent of time, indicating that the chain is pinned 
to one of its meta-stable states. In Fig. ^(b) the same 
chain is driven by a stronger force, F — 0.29. The dy- 
namics is characterized by jerky, non-uniform motion. In 
the time interval monitored, the chain moved very slowly 
over a distance of the order of its length. Increasing the 
force further, the trajectories become smooth again and 
are little affected by disorder. Figure [lj(c) shows an ex- 
ample with F — 0.4. The beads march forward with a 
finite velocity which is given by the slopes of the trajec- 
tories. By measuring the average slopes for different F's, 
we determine the velocity- force characteristics v(F). 

Numerical results for v{F) are shown in Fig. O for 
systems of size A^ = 1024. Depending on the value of F, 
the time averages were taken over the intervals of 10^ to 
10^ steps. Further disorder average was performed over 
10 to 50 different realizations of {Vj} and {un}- (For F 
far exceeding Fc, there is no need for large number of 
samples.) The data clearly indicated a sharp rise in v at 
a threshold force of Fc ~ 0.289. In the inset of Fig. [O, we 
plot V against the reduced driving force 6F = {F — Fc)/Fc 
on log-log scale. This yields the expected scaling form, 
V ^ (SF)'^ , with P « 0.25. For comparison, we also show 
the corresponding v — F characteristics for the case of 
uniform springs (with a„ = 10) in Figs. We again 



find critical depinning behavior, with Fc ~ 0.275 and 
/3 ~ 0.41. 
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FIG. 12. Trajectories of a system of A = 16 beads for 
three characteristic driving forces F. Size of time step is 
At = 0.01. 

Figs. ^ and |lj clearly show the difference between 
the uniform and random spring systems. In the case of 
uniform springs, our estimate for exponent f3 is compa- 
rable with thp result /3 = 0.45 ± 0.05 obtained by Myers 
and SethnaEj in their simulations of a one-dimensional 
automation model believed to be in the CDW univessal- 
ity class; it is also consistent with the FRG resultcl for 
the ID CDW (/? = 1/2) described in Sec. IV. A. For the 
random chain, our value of /3 is close to the one obtained 
in numerical studies of the dynamics of a_directed elas- 
tic string in 2— dimensional random mediacJ, which finds 
(3 = 0.24 ±0.1 for strong pinning and (3 — 0.34 ±0.1 for 
weak pinning. 
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FIG. 13. Average velocity i; as a function of the force F 
exerted on the random spring chain. Dotted line is the best 
fit to the form v = vo{F — Fc)^ where constants vo, Fa, and /3 
are fitting parameters. Inset shows a log-log plot of v vs the 
reduced force 5F — {F — Fc)/Fc\ straight line indicates the 
suggested scaling behavior. 
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FIG. 14. The v — F characteristics for a uniform spring 
chain. Dotted and straight lines are the best fit curves anal- 
ogous to the ones shown in Fig. 

Next, we characterize fluctuations in the bead con- 
figurations in the vicinity of the critical point. Posi- 
tional fluctuations imphcit in Figs. |l^ can be visual- 
ized more directly by plotting the displacement fleld 
Un{t) = rn{t) — Rn- In Fig. |l^, we show the typical 
response of the random chain with — 1024 beads to 
the applied force in the three regimes below/ at /above 
the depinning transition, all with the same realization 
of disorders. For F < Fc (Fig. |l^(a)), all beads stop 
to move after some transient motion characterized by 
small, local rearrangements. For F k, the mo- 
tion of the system is highly nonuniform as shown in 
Fig. |l^(b). (Here, the fines plotted are the displace- 
ment profiles for t = 0, 1000, . . . , 5000 with At = 0.01.) 
The profile advances in a very jerky fashion reminiscent 
of avalaaehes observed in models of sandpile and earth- 
quakegl^ca. Note that these avalanches have sizes com- 
parable to the system size, making the displacement pro- 
files much "rougher" than the ones shown in Fig. E5^(a). 
For F > Fc, the avalanches overlap each other and the 



displacement profile, cannot stop moving as they did in 
the previous casea£3. Fig. ^(c) shows the system evo- 
lution when the driving force is much larger than Fc- 
(The chain positions are plotted at a time interval of 
t — 1000 as in Fig. |l^(b)). Since this time interval much 
exceeds the avalanche overlap time, individual avalanches 
are not discernible at this scale, and the system advances 
smoothly witlp—a reduced roughness (governed by the 
KPZ equationN). 
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FIG. 15. Temporal evolution of the displacement profile 
for the random chain of A*' = 1024 beads at different driv- 
ing forces F. The lines plotted are displacement profiles at 
t = 0, 1000, . . . , 5000, with time step of At = 0.01. 

We now characterize the roughness of the displacement 
profile quantitatively close to the critical point. We moni- 
tor the disorder-averaged equal-time correlation function. 



(36) 



Systems with = 8 to 256 beads were examined right at 
the threshold forces F^- Simulations were run until the 
systems become "barely pinned", defined operationally 
as the point where v < 10"'^. Disorder averages were 
performed over 5000 samples for the smaller iV's and 100 
samples for the largest N. 

Numerical results for C(n, N) are shown in Fig. |l^ for 
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the random chain. Because of the SOS-Uke restriction 
we imposed on the local dynamics, C^^^(n, TV) can at 
most be linear in n as discussed in Ref. |2^. This upper 
bound is reached by all the curves shown in the figure. 
Thus C^/^{n, N) ~ n^, with % « 1. However, the data 
obviously contain additional dependence on the system 
size N and suggest the form C^/^(n, N) = n-N^'^^ , where 
Xfs is the finite-size exponent. To obtain this exponent. 




n 



FIG. 16. Correlation functions for random spring chains of 
size A'^. The straight line shows C{n,N) ~ |n|^. Inset shows 
the effective finite-size roughness exponent xfs(-^), defined 
by Eq. (^7|). Dotted line is the suggested asymptotic value of 
the exponent. 




,3 



n 



FIG. 17. Correlation functions for the uniform chains, 
plotted in the same way as Fig. |l^. 

we compute an effective exponent XFsi^) defined a 
xMN = T)^-L_±lo,,[^^$^^). (37) 

The result is shown in the inset of Fig. |l^. We see xfs (N) 
stabilize for N > 100, yielding xfs = 1-22 ± 0.01 for 
the random chain. The same calculations finds xfs = 
1.46 ± 0.01 for the uniform spring chain (Fig. [l7| ). 

The obtained value for xfs in the case of janiform chain 
is in good agreement with the FRG resulta which found 



Xfs = (4 — D -I- r]s)/2, with rjs = to one-loop or- 
der. It is also consistent with the numerical result of 
Xfs — 1-3 ± 0.3 found in Ref. The result for the 
random chain is also very close to the one reported in 
Ref. for the simulation of a driven elastic string with 
random-field or random-bond disorder (xfs — 1-25), and 
to the ones obtained from the simulations of related mod- 
els of interface depinning (xfs = 1-23 ± 0.02 in Ref. |9^ , 
and Xfs = 1.2 ± 0.1 in Ref. p6| ). Combining the re- 
sults on the two independent exponents /? and Xj we find 
strong evidence supporting our expectation that the crit- 
ical depinning dynamics of the driven random chain is in 
the same universality class as the driven elastic string in 
(1 -f- l)-dimensional random medium. This is the D = 1 
case of our more general conjecture concerning the crit- 
ical dynamics of the D-dimensional randomly-tethered 
elastic network. Numerical studies of the dynfuaics of 
the two-dimensional system is already underwayllj. 

V. BULK-MEDIATED ELASTICITY 

So far we have studied the properties of the randomly- 
tethered elastic network which is completely in contact 
with a disordered substrate. In many situations however, 
the elastic medium interacts with the substrate only at 
one of its surfaces. This is for example the case of friction 
between a thick piece of rubber and a piece of sandpa- 
per. Similar models have been used to describe aspects 
of "tectonic plate" movement along an earthquake fault 
zone (See Fig. |l|). Here we extend the model of Sec. IV 
to include bulk-mediated nonlocal elasticity, and study 
the driven dynamics of such systems. 

F 




FIG. 18. Picture illustrating two "tectonic plates" sliding 
past each other along a "fault zone". 

We shall focus on the (1 -|- l)-dimensional caricature 
of the problem depicted in Fig. |l^. We consider a 
two-dimensional random bead-spring network (of size 
Lx 'X Lz) in contact with a onc-dimensional disordered 
(and impenetrable) substrate V{x) located at z = (see 
Fig. |l^) . We wish to study the dynamics generated by a 
driving force applied to the upper [z — L^) surface. For 
simplicity, we assume there is a sufficiently large loading 
force —Fmz which keeps the elastic network in contact 
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with the substrate, thereby allowing us to suppress the 
z-degrees of freedom in the displacement vectors u. 




FIG. 19. A 2D bead-spring system of thickness Lz on 
disordered one-dimensional substrate. 

Let the a;-component of u(R„) at the upper and lower 
boundaries of the array be wq (n) and wi (n) respectively. 
Consider first the homogeneous elastic medium. The 
equations of motion for the w's are then 

Ho^dtwo = ^ {-G^^{n - n',0)woin\t) 



-G-\n~n',L,)wiin',t)} 
~V'{n ■ a + Woin, t)) 
J2{-G^\n-n',0)wi{n',t) 



(38) 



-G-\Ti-n',L^)wo{n',t)} + F (39) 



where /Zj's are the microscopic frictional coefficients, and 
the G{n,zys are the Green's functions whose Fourier 
transforms are given in Appendix B. We assume that 
the upper surface wi relaxes very quickly such that the 
nontrivial dynamics is dominated by the behavior of the 
lower surface wq, due to its contact with the substrate. 
Setting dtWi = and solving Eq. (^9|), we obtain a closed 
equation of motion for wq: 

fiodtwo = ^ J(n - n')wo{n',t) 

n' 

-V'{n-a + WQ{n,t))+F, (40) 

where the kernel J(n) is given in term of its Fourier trans- 
form 




E 



1 



Pn 



Pn = rnrL^ (41) 



n— — oo 

for Ifcl Lz :$> I 



-2-fL,k^ for \k\L^ < 1 ' 



(42) 



as detailed in the Appendix B. Thus, the familiar form 
of bulk- mediated elasticity, J{k) ^ — |fc| is recovered in 
the limit of large L^, while for small L^, the problem is 
effectively one-dimensional in the limit of small k. 

Randomness in the elastic medium can be readily in- 
corporated into the dynamics. Since the ground state 
configuration of the beads are unfrustrated, the effect of 



random spring lengths can be shifted away completely 
upon using the displacement field m(R„) defined with 
respect to the equilibrium position R„ — (2;„,z„) of the 
beads in the random system (without the external poten- 
tial V{x)). Using WQ{n) — uo{n) + Xn — n ■ a in Eq. (|4^), 
where uo(^) = u{xn, Zn — 0), we obtain finally the full 
equation of motion 



Mo ^^t"o = ^Ji'f^ ~ n')ua{n',t) 

n' 

-V {xn + uo{n,t)) + F. 



(43) 



The coarse-grained dynamics generated by the discrete 
system (^ ) with the kernel (|4^) can be derived as before. 
The resultant equation of motion is 



fj.Q^ dtUo{x,t) ^ / J{x - x')uq{x' ,t) 

J x' 

-V'{x)+Ui4x,uo)+F, (44) 

where the pinning force is /pin(a;,Mo) — V^{x)dpo{x — 
Uo,z = 0), with Spq{x,z) being density variation of the 
relaxed bead-spring system as described before, yielding 



[/pin(a::,Mo)/pin(0, 0)] w A6{x)S{uo), 



(45) 



for a variety of randomly-tethered networks. 

Eq. ( p4[ ) with the disorder correlator (45) is very simi- 
lar to the equation of motion of a contact line which con- 
trols the spreading of a non- wetting liquid on a disordeced 
solid surfacel£3. As shown by Joanny and de GennesEfl, 
a contact line is governed by nonlocal elasticity of the 
form (^2|), reflecting the energetics of distorting the un- 
derlying liquid/gas interface. The driven dynamics of tlie| 
contact line has been investigated by Ertas and KardarO 
using the FRG method. A depinning transition similar to 
that of the driven KM was found. Upon generalizing the 
substrate to i)-dimensions, an expansion about the up- 
per critical dimension D = 2 yields the exponent values 
X = e/3 and z = 1 — 2e/9 to first order in e = 2 — D. Other 
exponents can be obtained from the exponent relations 
jy = 1/(1 - x), P ={z - x)/(l - X), andK = D-l + x- 
Thus the contact line problem (e = 1) is described by 
the exponents w 3/2, /3 « 7/9, and x ~ 1/3- Our 
problem of course differs from that of the contact line 
again by the long-range correlations in higher moments 
of the pinning force f{x,uo). However, as in the case of 
the random spring with local elasticity, we do not expect 
these higher moments to change the universality class of 
the depinning dynamics. 

To test this conjecture, we performed numerical inte- 
gration of the discrete equation of motion (p3[), with the 
effective kernel J precomputed using (^H), with — 128 
and Lz = 4096. (Similar results were obtained when 
we directly used the kernel J — — 27|fc|.) The parame- 
ters describing the springs and the random potential V 
were the same as the ones used in Sec. IV. B. The data 
for the velocity-force characteristic (shown in Fig. [20|) 
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were obtained by averaging over 15 independent sam- 
ples, while the spatial correlation function of uo{n) at 
F = Fc (Fig. |2l|) was averaged over 100 samples. From 
the numerical data, we find the exponents (3 ~ 0.66 and 
X ~ 0.25, which are consistent with the approximate 
FRG results for the contact line. Thus the correspon- 
dence between the critical dynamics of the random spring 
chain and that of the directed path appears to hold even 
with nonlocal elasticityll2E . 
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FIG. 20. The velocity-force characteristics in the vicinity 
of the depinning threshold, for the 1-1-1 dimensional system 
of size 128 x 4096. Log-log plot is shown in the inset. The 
slope of the straight line gives (3 = 0.66. 
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FIG. 21. Roughness of the displacement field at the 
threshold Fc = 0.1910. The full line gives x = 0.25. The 
statistical uncertainties are smaller than symbol size. 

We have so far discussed only the critical dynamics 
of the elastic medium driven at a constant force. An- 
other way the system may be driven is by a constant 
velocity imposed at the upper boundary. Such [Situa- 
tions are of interest to the study of earthquakesEJ, as 
they model tectonic plate motiop-along a fault. The well- 
known Burridge-Knopoff modeS is of this class. 

To consider the effect of a constant-w drive, we re- 
turn first to the linear bead-spring chain described in 
Sec. IV. B. We drive the chain by connecting each bead 
n to a reference point i?* via a weak loading spring of 
spring constant a. The spacing of the reference points 
are chosen such that when the beads are in the relaxed 



state (i.e., without the external potential V), the loading 
springs are unstretched. Thus, i?* — i?„ = X]m=i 
where the lengths of the springs connecting the 

beads. The reference points are then set to motion, with 
i?*(t) = Rn + vt. This motion forces the bead-spring 
chain to move, at the same velocity u, via the action of 
the loading springs (see Fig. 



R'L (0 




Unit) 



FIG. 22. A one-dimensional random spring chain pulled 
by a constant velocity drive via a set of loading springs of 
spring constant a. The mean lag distance is u. 

The equation of motion for the beads is 

dtTn = l[rn+i - 2r„ -1- r„_i - (a„+i - a„)] 
-a[r„-i?;(t)]-l/'(r„), 

where r„(i) is the position of the bead n. In term of the 
displacement field Un{t) = rn{t) — R^{t) + u, {u is the 
mean lag distance as shown in Fig. the equation of 
motion becomes 

-v + au-V'{Rn + vt + Un{t) -u). (46) 



It is instructive to compare Eq. ( 
uated in the moving frame r„ (t) 



|) with Eq. (PI) eval- 
Rn + vt + u^Jt): 



dtUn = 7 [""n+l - 2u„ + Un-l] 

-V-V'{Rn+Vt + Un{t)) 



F. 



(47) 



We find the systems described by ( ^ ) and ( ^ ) to be 
statistically the same, up to the damping term — q;u„, if 
we identify au in ( p6| ) with the driving force F. (The 
random forces V'{x) and V'{x — u) clearly have the same 
statistics.) The knowledge that Eq. ( p7| ) exhibits a de- 
pinning transition with v (x {F — Fc)~and a diverging 
correlation length ^ ~ {F — Fc)~'^ immediately lead us 
to conclude that large fluctuations also occur in the sys- 
tem with constant-w drive as u — > 0, with ^ ~ v~'^^^ 
and u — > Fc/a. However, true critical behavior is pre- 
vented by the linear damping term —aun in (Hw. The 

1 /'^ — 

correlation length saturates at ^* ~ (7/a) • 

We now turn to the system depicted in Fig. |l^, with 
the upper (z = Lz) boundary set to constant velocity. 
Mathematically, the situation is described by Eq. (|3^), 
with zi'i(n,t) = vt. The explicit form of the dynamics 
is readily obtained from the knowledge of (see Ap- 
pendix B); we find the following equation of motion. 
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dtU{){n^ t) — J{n — n!) UQ{n\ t) 

n' 

-aL{uQ-vt) -V'{xn+uo{n,t)) (48) 

where the kernel J{n) is approximately the same as that 
describing the constant-/ dynamics in Eq. (^3|), and the 
damping coefficient is ~ y /(nLz). In the moving 
frame uq ^ uq + vt — u, Eq. ( [48| ) has the same form as 
( p6| ) except for the nonlocal elasticity coupling J. Thus, 
we expect similar "near-critical" behavior for the sys- 
tem (48) as V —>■ 0, with a diverging correlation length 
^ ~ v'"/'^ which saturates at ^* « (7/aL)^/'^. The value 
of the exponent a depends on which regime of J the sys- 
tem is in at the scale see Eq. (|4^). We have a — 2 li 
Lz/£,* 1 and cr = 1 if L^/i* » 1. Given the expres- 
sion for aL^Lz), we find cr = 1. Thus, ^* ~ Lz becomes 
very large and the system is very nearly critical as u ^ 0. 
(Note that the lag distance u also becomes very large.) 
Zero-temperature driven dynamics with w — > studied 
here is an example of "extremal dynamics" studied ex- 
tensiyely in the context of self-organized critical phenom- 
enaEJ. The relation between self-organized criticality and 
dynamic critical, plieuoHiena has been discussed in vari- 
ous context^eHEle. Within our model spring-bead 
system, we see that the equivalence between these two 
classes of phenomena can be explicitly established. 

We finish this section with a discussion of the {D + 1)- 
dimensional generalization of the system depicted in 
Fig. ^ Our analysis lead us to anticipate strongly 
that such a system is equivalent to the _D-component, 
_D-dimensional driven manifold in (D + £))-dimensional 
random media, with a nonlocal elasticity kernel (^2|). In 
particular, the D ~ 2 case is analogous to the sliding 
slabs of elastic media (Fig. |l8|) . Since _D = 2 is the crit- 
ical dimension for problems with bulk-mediated elastic- 
ity, the critical[£xponents at the depinning transition are 
known exactlytl, e.g., x = O(log), v = \, and (3 = 1. 
From these and the exponent relations given above, one 
finds the exponent describing the avalanches at the onset 
of motion to be k = 1. Recently, the avalanche distcibu- 
tion of a related earthquake model was investigatedEiJ. In 
the study of earthquakes, one monitors the distribution of 
the "moment" M = J d^x Au(x), where Au(x) is the to- 
tal displacement at location x during an avalanche. Since 
X = in D = 2, M ^ P where £ is the linear size of the 
avalanche, assuming that the avalanche clusters are com-, 
pact up to logarithmic correction. From (|2^), it foUowsE^ 
that the probability Pr(Af ) of finding an earthquake with 
moment exceeding M is Pr(Af) ^ M~^p{M/(^) at the 
critical point, with B = k/2 = 1/2. The model of Ref. |5l| 
uses constant-w drive with v ^ 0. Thus, 1^ = ^* ~ Lz- 
The numerical result on the moment-distribution is con- 
sistent with the RM analogy, as has been pointed out in 
Ref. Fl|- 



VI. SUMMARY 

We have presented a detailed study of the static and 
dynamic properties of an elastic medium interacting with 
a disordered substrate. The behavior depends crucially 
on the equilibrium density distribution of the elastic 
medium in the relaxed state. The interaction of a per- 
fectly homogeneous medium with the substrate belongs 
to the CDW universality class as has already been dis- 
cussed in the literature. The somewhat surprising re- 
sult of this study is that a slight amount of quenched- 
in inhomogeneities of the elastic medium, even with 
only interstitials/ vacancies or phonons, is sufficient to 
change the universality class. Instead of the CDW, a D- 
dimensional inhomogeneous medium on Z?-dimensional 
disordered substrate belongs to the universality class of 
a I?-dimensional homogeneous manifold embedded in an 
effective {D + £))-dimcnsional random medium. This 
is a consequence of the fact that quenched-in density 
variation breaks the translational symmetry of the elas- 
tic medium, such that the dense/dilute parts of the 
medium preferentially stick to the attractive/repulsive 
parts of the substrate. We verified this numerically for a 
one-dimensional random bead-spring system: The finite- 
temperature static behaviors of the random chain are 
found to be indistinguishable from those of the directed 
path in (1 -I- l)-dimensional random medium. The zero- 
temperature driven dynamics exhibits a depinning tran- 
sition, the critical properties of which are also indistin- 
guishable from those of the driven elastic string in 1-1-1 di- 
mensional random medium. The equivalence is found to 
hold also for elastic media with nonlocal (bulk-mediated) 
elasticity, making our model system and results relevant 
to tribology problems. Finally, a slow- velocity drive is 
shown to be equivalent to a constant-force drive close to 
the depinning transition, demonstrating explicitly that 
self-organized critical phenomenon obtained from certain 
extremal dynamics may be viewed as a dynamic critical 
phenomenon. 

Although we formulated the interaction of the elas- 
tic medium with the substrate in terms of the energetics 
of density variations, the underlying physics governing 
the interaction is much more general. The situation be- 
ing studied here is really one of pattern matching, i.e., 
matching of quenched-in fluctuations between two differ- 
ent elastic media. This could occur as well in the form 
of roughness matching, say, the matching of two rough 
surfaces in contact, or sequence matching, as in the hy- 
bridization of two heterogeneous DNA sequences. By 
viewing interacting random systems as effective homoge- 
neous systems embedded in higher spatial dimension with 
external randomness, it is possible to simplify and resolve 
a large class of interesting problems. These include for 
example ttie reptation of heteropolymers in disordered 
gel matrixa, where "pattern matching" of the polymer 
composition with its reptation tube leads to anomalously 
slow dynamics, and the preferential adsorption of het- 
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eropolyraecs on surfaces coated with specific chemical 



patterns 



We hope that this work wiU stimulate fur- 



ther progress in understanding the physics of interacting 
random systems. 
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Note added: After the submission of this manuscrict, we 
became aware of recent works by K.V. SamokhinE£3, who 
considered theoretically the case of an amorphous mani- 
fold (with local elastcity) on random substrate. Conjec- 
ture on the perturbative irrelevancy of long-range corre- 
lated effective random potential was asserted based on a 
replica calculation. 



APPENDIX A: COARSE-GRAINED DYNAMICS 

In this appendix we derive the coarse grained equa- 
tion of motion (|8|) (with the definitions (|2|) - (^), 
starting from the discrete model of the inhomogeneous 
elastic network driven by a constant force. For simplic- 
ity, we shall describe only the case D = 2, and further- 
more exclude all quenchcd-in topological defects includ- 
ing interstitials and vacancies. Inhomogeneities in the 
medium are described by small deviation w„j^_„2 of the 
equilibrium positions of the beads from a perfect lattice. 



Tin 



niai -|- n2Si2: where the a^'s are the lattice vec- 



tors. For |w| <C |a|, we can use n = (711,712) to label the 
beads. 

We start with the coarse-grained density field, 



where 



Pir,t) = ^ Jj^r' pir + r') 



p(r)=^5(r-r„(t)) 



(Al) 



(A2) 



is the microscopic density field, rn{t) gives the actual 
position of the bead n, and A is the coarse-graining scale, 
of the order of several |a|'s. In Sec. II, we showed that in 
terms of the displacement field u(r). 



p(r) « p(l - V • u) + 6po{r - u(r)). 



(A3) 



where p = l/|ai x a2| is the average bead density, and 
i5po(x) describes the equilibrium density variations. For 
the specific 2D model considered here, 



Spo{x) - p^cos[Kj • (x + w(x))], 



(A4) 



with Kj being the reciprocal lattice vectors. Since Spo{'x.) 
fluctuates predominantly at the scale |a|'s, it does not 
survive the coarse graining, and we have 



p«p(l- V-u) 



(A5) 



to leading order in (V • u). 

The equation of motion for u can now be obtained 
from the evolution of p, which is given by the continuity 
equation, 

9tp(r,t)-HV-J(r,i) =0, (A6) 
where j is the coarse grained "current" , given by 

J(r,i) = ^^dVj(r-|-r',i), (A7) 

and 

j(r,t) =^9ir„5(r-r„(t)). (A8) 

n 

This choice of the current automatically satisfies the con- 
tinuity equatio n, a s c an be verified directly by inserting 
the definitions (^)-(^^d (|A^-(^) into Eq. (^). 
Using (A5) for p in Eq. (A6), we easily find the form of 
the equation of motion for u: 



dtu{r,t) = j(r,t)/p. 



(A9) 



Thus our task is to obtain the coarse-grained current j 
starting from the discrete equation of motion, 

fJ^O^dtTn = ^ 7(l"n+a - fn + Wn+a - W„) + f (r„), (AlO) 



where the sum is over the nearest neighbors, and 

f (r) = - VF(r) -t- F (All) 

is the external force. 

Consider first the "one-particle" equation of motion, 
PQ^dfYn = f(r„(t)). Then from the definition (A7), we 
have 



j(r,i) = fiap{r,t)i{r) 

= /io(l- V-u)f(r) + Aiofi(r,u(<)), 



where 



fi(r,u) = Spo{r - u) • f(r)/p. 



Upon coarse graining of Eq. ( A12 ), we find 

j = /io(l - V • u)f(r) +/iofpin(r, u(<)), 
where 



(A12) 
(A13) 
(AM) 
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f (r) = A-^ [ (fr' f (r + r') = - VF< (r) + F (A15) 
Ja 

is the coarse-grained force, V<(r) being the slowly vary- 
ing part of the substrate potential, and 



fpin(r,u) = -Vl/>(r) <5po(r-u) 



(A16) 



(r) denoting the Fourier modes close to the reciprocal 
lattice vector K^'s. 

Inclusion of bead-bead coupling as specified by 
Eq. ( AlO| ) only affect the term f. For a statistically 
isotropic array (i.e., a triangular lattice), we find 



f =(cn(r)-C66(r))V(V-u) 

+C66(r)V2u- Vy<(r) +F 



(A17) 



where the elastic moduli cii,C66 ~ 70^. (The r- 
dependences of the c's come from V-w(r) and V x w(r).) 
The equation of motion for the displacement field u is fi- 
nally 



^i^^dtU = (1 - V • u) f + fpin(r, u) 



(A18) 



with f and fpi„ given by Eqs. (A17) and ( |A16| ) respec- 
tively. 



APPENDIX B: BULK-MEDIATED ELASTICITY 

In this appendix, we derive the nonlocal equation de- 
scribing the motion of the beads at the z = boundary of 
an elastic medium placed in a; — z half-plane (see Fig. 19). 
Our strategy is to obtain first the effective Hamiltonian 
describing the equilibrium fluctuation of the beads at the 
z = Q boundary, subject to various drive conditions ap- 
plied to the opposing (2 = L) boundary. The effective 
equation of motion is then obtained by applying gradient 
descent dynamics. 

We consider here the homogeneous elastic medium 
which is described by the Hamiltonian (^. To simplify 
the description, we assume that a sufficiently large nor- 
mal force is applied such that displacement in the trans- 
verse {z) direction is negligible. The Hamiltonian de- 
scribing the residual displacement fluctuation w{x, z) in 
the x-direction can then be written as 



]^ J dxdz^{Vw{x,z)f . (Bl) 



mo [u 



The probability to find a system in the state with dis- 
placement WqIx) at the z = boundary and z«i(x) at 
the z = Lz boundary is 



P[wt),wi\ oc j Vw 5 {w{x,0) — wq{x)) 

xS{wix,L,)-wi{x)) e-^^°^'"\ (B2) 
Integrating out the w-field, we obtain 



P{wq, wi} oc e 



(B3) 



where 



- -G-i(fc,o) iz5)o(fc)r + R(fc)r 



-G-\k,Lz)wo{k)wi{-k)} , (B4) 
with Wi{k) being the Fourier transform of Wi{x), 

C(fc,0) 



G-\k,0) = 
G-\k,Lz) = 



C2(fc,0)-C2(fc,L,) 

C{k,Lz) 
C^{k,Q)-C^{k,Lz) 



(B5) 
(B6) 



and 
C(/c,z) 



E°° cos(pz) 



2^Lz ^ fc2 + p 

P— — OQ 

-\k\z 



o,±^,. 



27l*:| 



for \k\L, > 1 



27L, 



i^^ (1 + fc2i2/7r) for |fc|i^<l 



(B7) 



The equations of motion for and w\ are now straight- 
forwardly obtained from gradient descent of PH-cs- For 
a constant force drive applied to the z = boundary, 
we have 

fio'dtWoik,t) = -G-\k,0)ido{k,t) 

-G-^{k,Lz)wi{k,t) (B8) 
fi^^dtWi{k,t) = -G^^{k,Lz)wo{k,t) 

-G-\k,0)wiik,t)+FS{k) (B9) 

where fii are the different microscopic frictional coef- 
ficients at the boundaries. Assuming that the upper 
boundary wi relaxes very quickly such that the nontriv- 
ial dynamics is dominated by the lower boundary (due 
to its contact with the substrate), we can set dtwi — 
and obtain an effective equation for uq itself. It reads in 
Fourier space 



with 



fiodtWo{k, t) = J{k)wo{k, t) + F 5{k) 



J{k) ^ -C~i(fc,0) 

-2j\k\ for |fc|i^>l 
-2-fL,P for \k\Lz < 1 



(BIO) 



(Bll) 



Thus, in the limit of large Lz, the kernel J takes on the 
form J{k) = —2j\k\ familiar for bulk-mediated elastic- 
ity. And in the opposite limit, the coupling becomes 
local again as the system reverts back to being one- 
dimensional. 

Equation of motion for constant w-drive at the z = Lz 
boundary is obtained simply by setting wi (fc, t) = vtS{k) 
in Eq. (|B8[). It reads in Fourier space 
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dtWo{k,t) = -G-^{k,0)wo{k,t) 

+ G-\k,L,)vot6ik). (B12) 

Using the result (p^) in Eqs. and @, we find 



for |fc|L^>l . . 



^ '^'''^^ ~ {«! +7Li^.fc^ for |fc|L, « 1 
G-\0,L,)=aL 



(B14) 



witli aL ~ lliT^Lz) and 7^ oc 7. Tims the equation of 
motion for the constant-w drive is 

dtwoik, t) « -aL {wo{k, t) - vt5{k)) 

+J{k)wo{k,t), (B15) 



with J(fc) given approximately by Eq. (Bll) 
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